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CO . ABSTRACT 



Context. The low-mass star formation evolutionary sequence is relatively well-defined both from observations and theoretical consid- 
erations. The first hydrostatic core is the first protostellar equilibrium object that is formed during the star formation process. 
Aims. Using state-of-the-art radiation-magneto-hydrodynamic 3D adaptive mesh refinement calculations, we aim to provide predic- 
^0 ' tions for the dust continuum emission from first hydrostatic cores. 

Methods. We investigated the collapse and the fragmentation of magnetized 1 M Q prestellar dense cores and the formation and evo- 
lution of first hydrostatic cores using the RAMSES code. We used three different magnetization levels for the initial conditions, which 
cover a wide variety of early evolutionary morphology, e.g., the formation of a disk or a pseudo-disk, outflow launching, and fragmen- 
tation. We post-processed the dynamical calculations using the 3D radiative transfer code RADMC-3D. We computed spectral energy 
t-H ' distributions and usual evolutionary stage indicators such as bolometric luminosity and temperature. 

| Results. We find that the first hydrostatic core lifetimes depend strongly on the initial magnetization level of the parent dense core. We 

derive, for the first time, spectral energy distribution evolutionary sequences from high-resolution radiation-magneto-hydrodynamic 
calculations. We show that under certain conditions, first hydrostatic cores can be identified from dust continuum emission at 24 fim. 
and 70 fim. We also show that single spectral energy distributions cannot help in distinguishing between the formation scenarios of 
the first hydrostatic core, i.e., between the magnetized and non-magnetized models. 

Conclusions. Spectral energy distributions are a first useful and direct way to target first hydrostatic core candidates but high- 
resolution interferometry is definitively needed to determine the evolutionary stage of the observed sources. 

in ■ 

\Q ■ Key words. Magnetohydrodynamics (MHD), radiative transfer - Methods: numerical- Stars: low mass, formation, protostars 
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' 1. Introduction tions may n ot always reflect the true evolutionary stage of an 

CN '■ object (e.g., iDunham et aD2010l) . 

It is well established that low-mass stars form from the col- 

lapse of dense prestellar cores. The star formation process can The protostellar collapse phase is divided into two phases . 

• H. be dlvlded both observationally and theoretically into an evo- Fim the prestellar dense core ( grav itationally bound, no proto- 

V . lutionary sequence. From observations, protostars are classified star> and no QUtflow) collapses isot hermally and is able to freely 

^ ; according to the parametrized properties of the observed spectral radiate itfj gravitational energy int0 spacCi until it becomes dense 

Cd , energy distribution (SEP) and appearance. For instance, a Class ^ „ 10 -i 3 cm - 3) and opaque enough for that radiation t0 

source d Andre et al. || 1993|) , which corresponds to the phase be trapped At this ^ the gafj beg ins t0 heat up and the first 

where the mass of the protostar is less than half the mass of hydrostatic core (FHSC) is formed jLarsonlll969h . The FHSC 

the surrounding envelope, is defined according to its bolometric is optically thick and its temperature increases adiabatically un- 

temperature T bol d Chen et alj|1995|) or to the ratio of bolomet- til k reaches „ 20QQ K where Ro dissociation starts . This en - 

nc to submilhmeter luminosities L bol /L smm ( |Andre et al. || 1993|) . dothermic reaction allows the gas t0 underg0 a second collapse 

At later stages, the Class I, Class II, and Class III sources can phase unt0 stellar densities are reac hed, i.e., the formation of 

be classified according to the slope a IR = dtog(^)/dlog(^) the second hydrostatic core (SHSC), the protostar that corre- 

between 2.2 /zm and 10 - 25 fim ( |Lada| | 1987|) . These classi- sponds to the prot ostellar Class evolutionary stage. The FHSC 

fications attempt to replicate the different evolutionary stages ifj thufj a tmmition phase from pres tellar cores to Class objects 

of a protostar. Observational ambiguities (e.g., projection ef- and has not been identified conclusively because its obser- 

fects, episodic accretion), however, imply that these classifies yational appearance is unclear The pred icted timescales of the 

FHSC phase range from a few hundred years in spherical models 
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low-m ass non-magnetized rotating dense cores (ITomida et al.l 
I2010ah . 

In recent years, significant progress has been achieved in 
both theory and observation thanks to increasing computing 
power and new observational capabilities provided by Spitzer 
and Herschel. For instance, th e large Spitzer "C ore to Disks" 
(c2d) Legacy Program survey dEvans et alJl2003h helped to in- 
crease the number of known protostars up to the point where 
first detections of FHSC can now be expected, given the short 
durat ion of the FHSC phase dEnoch et al.ll2010t iDunham et alJ 
1201 lb . Spitzer also helped to identify a new class of low- 
lumi nosity objects: the Very L ow Luminosity Objects (VeLLOs, 
e.g., Idi Francesco et al.ll2007l) . These are embedded in a dense 
envelope and have low intrinsic luminosity, L < 0.1 L Q . It is still 
a matter of debate whether VeLLOs are dense cores undergo- 
ing first collapse or lo w-mass pro t ostars with a low-accretion 
rate. In a recent study, iKim et aL I (l20Tlh showed that at least 
in one case (the Bok globule CB 130-1), the latter is the most 
likely explanation. Prestellar dense core populations have been 
identified over the pas t twenty years using space- and ground- 
based telescopes (e.g., Ward-Thompson et al.ll 1 994t iMotte et al] 
119981: iJohnstone et all 120001 iNutter & Ward-Thompsonll2007l) . 
More recently, Herschel observations have i dentified hundreds 
of prestellar cores and Class sources (e.g. , I Andre et"a"Dl2010t 



iBontemps et al.ll2010t Henning et ail 2010h . Combining differ 
ent instrument data. lLaunhardt et al.l d201 0ft presented the results 
of a comprehensive infrared, submillimeter, and millimeter con- 
tinuum emission study of isolated low-mass star-forming cores, 
and identified nine Class sources. Finally, a few FHSC can- 
didate detections have been reported in recent years thanks to 
the increasing spectral and spatial scale coverage ( Belloche et alJ 
2006t IChen et al.ll2010t lEnoch et al]|20Tot IDunham et al.ll201 It 



Pineda et alJl20 lit" ?). Their confirmation remains controversial, 
however, because the expected observable signatures of FHSC 
are still uncertain. 

In parallel, enormous work has been put into developing nu- 
merical hydrodynamical models that include magnetic fields and 



radiat i ve transfer (e.g., Price & Batd l2009t ICommercon et alJ 



1201 Ot ITomida et ai.1l20l6bF but relatively little has been done 
in producing synthetic observations that are directly comparable 
with actual observations. Most of the observables used to con- 
front theory and observation remain statistical in nature, such 
as the initial mass function (IMF) or the binary distribution, 
whereas observations report mostly dust continuum and molec- 
ular line emission data. Therefore, it is vital that the interpreta- 
tion of the data is sound, necessitating the use of non-symmetric 
models incorporating complex physics to derive predictions of 

dust continuum and molecular line emission obse rvations. 

Fo r the early stages of star formation, iBoss & Yorke 
(1 1995 1 ) . iMasunaga et al] (1 19981). lOmukail (120071). lYamada et al. 
d2009l). ITomida et all d2010al) . iTomisaka & Tomidal d201 lb . and 
ISaigo & Tomisakal ( 20111) made FHSC observational predic- 
tions of either dust emi ssion or molecular line emission. 
lYoung & Evansl d2005l) and lDunham et alJ d2010h also presented 
evolutionary signatures of the formation of low-mass stars, but 
used a crude approximation of the FHSC and the evolution of 
the envelope during collapse. To date, there has been no sys- 
tematic study of observational predictions for FHSCs that ex- 
plore physical conditions, and in particular on the magnetization 
level. Magnetic fields have indeed been f ound to control the first 
collapse and fragmentation pha ses (e.g., iHennebelle & Tevssierl 
120081: ICommercon et al.|l2010h and strong magnetization seems 
to be a favored sce nario for reproduci ng observational data of 
Class multiplicity dMaurvet alJl2010l) . The ability of a disk to 



fragment depends ion the magnetization level, leading to poten- 
tial changes in the SED. 

In this paper, we derive SED evolutionary sequences and 
classical observational indicators such as Lboi and Tboi from 
state-of-the-art radiation-magneto-hydrodynamic (RMHD) 3D 
calculations of collapsing 1 M dense cores with different initial 
magnetic field strengths. We address the questions of whether 
or not SEDs can help in distinguishing physical conditions (e.g., 
magnetic field strengths) and how SEDs can help in targeting 
FHSC candidates. This paper is the first of a series and is asso- 
ciated with a companion paper in which we produce synthetic 
ALMA dust continuum emission maps (Commercon et al. in 
prep, hereafter Paper II). 

The paper is organized as follows: in Sect. 2, we present the 
numerical codes we used and our post-processing methodology. 
The RMHD calculations are qualitatively presented in Sect. 3. In 
Sect. 4, the RMHD calculations are post-processed to produce 
synthetic SEDs. In Sect 5., we discuss our results and potential 
observational tests. Sect. 6 concludes our paper. 

2. Method 

2.1. The RAMSES code 

For the dynamical RMHD calculations, w e used the adap - 
tive mesh refinement (AMR) code RAMSES dTevssieri |2002|) 
which integrates the equa tions of ideal MHD dFrom ang et al. 
I2006b iTevssier et al.l I2006I) and the equations of radiation hy- 
drod ynamics under the gray fl ux-limited diffusion approxima- 
tion dCommercon et al.ll201l3) . The RAMSES code is well de- 
signed for collapse calculations and has been extensively used in 



low-mass star formation st udies ( e.g.. IHennebelle & Fromang 
2008t ICommercon et al.l l2008b IHennebelle & Ciardil 120091 
Commercon et a 1J 120101) . It was als o recently used in the 



high-mass star formatio n framework dHennebelle et al.l 1201 It 
ICommercon et alJl201 lbl) . 



2.2. The RADMC-3D code 

To post-process the RMHD calculations, we used the 3D ra- 
diative transfer code RADMC-3lQ, developed by C. Dullemond. 
RADMC-3D can work with 3D AMR grids such as the RAMSES 
grid. RADMC-3D has the capability to perform dust continuum 
calculations with anisotropic scattering and gas/molecular line 
transfer calculations (LTE or non-LTE). The dust temperature, 
emission, and absorption are computed using a thermal Monte 
Carlo method. RADMC-3D is able to compute images and spectra 
from dust continuum radiative transfer. For non-LTE line radia- 
tive transfer calculations, RADMC-3D can estimate the level pop- 
ulations using approxima te methods, such as the Large Velocity 
Gradient (LVG) method dShettv et al.l 1201 ll) . In this paper, we 
will only use the dust continuum features of RADMC-3D and ne- 
glect scattering. 



2.3. Interface 



First, we assumed that the dust and the gas are well cou- 
pled at high dens ities (p > 10~ 19 g cm -3 or n > 3 x 10 4 cm" 3, 
iGalli et all 120021) 



to estimate the dust temperature. Second, 
we assumed that the correct temperature is the gas tempera- 
ture computed in the dynamical RMHD calculations using the 
gray FLD approximation method. This second assumption is 



http://www.ita.uni-heidelberg.de/ dullemond/software/radmc-3d/ 



B. Commercon et al.: Synthetic observations of first hydrostatic cores in collapsing low-mass dense cores. 



3 




lier than the proto star+disk stage wher e SEDs peak at shorter 
wavelengths (e.g., iDunham et al.ll2010h 

3. Radiation-magneto-hydrodynamics models 

3.1. Initial conditions 



10.0 100.0 1000.0 10000.0 
Wavelength [/Lim] 

Fig. 1. Top: Scatte ring and absorption coefficient of 
ISemenov et al. I J200l as a function of the wavelength for 
the homogeneous spheres model (with Fe/Fe+Mg=0.3 "nor- 
mal" silicate composition). Bottom: Corresponding ratio 
between the scattering to absorption coefficients as a function of 
wavelength. 



likely co rrect for FHSC form ation and evolution, as shown re- 
cently bv lVavtet et alJ (l2012h for ID spherically symmetric gray 
and multi-frequency dense core collapse calculations. To obtain 
more accurate temperatures, we could have used the thermal 
Monte-Carlo capabilities of RADMC-3D. During the FHSC stage, 
however, the only sources of photons are from the compressional 
work and the accretion luminosity, which are not trivial to extract 
from the RMHD calculations. 

Thanks to the AMR grid capabilities of RADMC-3D, we di- 
rectly loaded the AMR structure of RAMSES within RADMC-3D. 
We thus did not lose any information between the RMHD and 
the post-processing radiative transfer calculations. We show in 
Appendix [A] the importance of using the AMR grid instead of 
remapping the data onto a uniform grid. 

We used the low-temperature opacities of ISemenov et alJ 
(I2003h for the homogeneous sphere model (with Fe/Fe+Mg=0.3 
"normal" silicate composition). For the RMHD calculations, we 
used a parametrization of the gray Rosseland and Planck means 
as functions of the gas temperature and density. For the post- 
processing radiative transfer calculations, we used the corre- 
sponding frequency-dependent opacities. Figure[T]shows the ab- 
sorption and scattering opacities we used, as well as the ratio 
of scattering to absorption coefficients. The scattering coeffi- 
cient is always less than the absorption coefficient, while their 
ratio is greater than 10% for visible to long near-infrared wave- 
lengths (i.e., 0.2 pm < A < 8 /mi). Since we neglected scat- 
tering, we may have overestimated the amount of radiation that 
escapes at shorter wavelengths. Neglecting scattering remains a 
good approximation, however, since in this study we focus on 
the internal luminosity of collapsing dense cores at a stage ear- 



We used similar initial conditions as in ICommercon et al.l 
d2010h . We considered 1 M Q uniform density and temperature 
sphere in solid body rotation about the z-axis. The initial radius 
of the sphere is R = 4.93 x 10 16 cm (~ 3300 AU). The initial an- 
gular velocity is set using a ratio of the rotational to gravitational 
energies B = 0.045. To favor fragmentation, we introduced an 
m-2 azimuthal density perturbation with an amplitude of 10%. 
The corresponding free-fall time and orbital time are tg ~ 33 
kyr and f 01 -b ~ 1.85 x 10 2 kyr, respectively. The initial density is 
po = 3.97xl0~ 18 g ctrT 3 and the initial temperature is To — 1 1 K. 
The ratio of specific heats y was set to 5/3. In this study, we did 
not include any additional microphysics such as H2 dissociation. 
Therefore, our models are only able to follow the dynamics up 
to a temperature of 2000 K, at which the second collapse starts. 

The magnetic field is initially uniform and parallel to the 
rotation axis. The strength of the magnetic field is expressed 
in terms of the mass-to-flux to critical mass-to-flux ratio p = 
(Mo/<t>)/(Mo/<t>) c . In this study, we explore three magnetization 
degrees: p = 200, which corresponds to a quasi-hydrodynamical 
model (model MU200), p — 10 for an intermediate magnetic 
field strength (model MU10), and p = 2 for a strong magnetic 
field case (model MU2). 

We used a coarse grid of 64 3 cells with ten levels of refine- 
ment, which corresponds to a maximum effective resolution of 
0.2 AU . Calculations were perfor med using the HLLD Riemann 
sol ver dMivoshi & Kusan otf2005h Following up on former stud- 
ies dCommercon et al.ll2.Q08l 20101) . we imposed at least 12 cells 
per Jeans length as a grid refinement criterion. 

3.2. Qualitative description 

In the following, we describe only qualitatively the results 
of the three RMHD calculations since the focus of this paper 
is the observational characteristics of the models and not their 
physical quantitative description. We synchronized the calcu- 
lations at time to at which the maximum level of refinement 
is reached, which corresponds roughly to the formation of the 
FHSC. Column density and temperature maps for the three cal- 
culations are shown in Fig. [2] All fragments that are formed in 
the calculations are FHSCs. 



3.2.1 . The strongly magnetized case - p = 2 

In the strongly magnetized case, fragmentation is totally sup- 
pressed and no rotationally supported disk is formed b ecause of 
the efficient magnetic braking (e.g., iJoos et aDl2012h . At time 
t + 0.78 kyr (with t = 50.4 kyr), there are 1.66 x 10 6 cells in 
the whole computational domain. An outflow is launched along 
the rotational axis with velocity ~ 3.7 km s _1 and it extends up 
to ~ 550 AU (not seen in the zoom in Fig. |2|. Only one object 
is formed, surrounded by a thick pseudo-disk region. The mass 
of the FHSC is ~ 0.06 M and its mean temperature is ~ 375 
K (T m . dx ~ 630 K). The heated region around the central FHSC 
is almost spheroidal, which cannot be obtained in models using 
a barotropic approximation. In the latter case, the warm region 
would correspond to the dense region (the temperature is set b y 
the density), i.e., the pseudo-disk (see lCommercon et alJl20l6T) . 
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MU2 , t=t +0.78 kyrs MU10 , t=t +1.81 kyrs MU200 , t=t +3.26 kyrs i og (N) 




-200 200 -200 200 -200 200 

x [AU] x [AU] x [AU] 



Fig. 2. Column density and temperature maps in the xy-plane (two upper rows) and xz-plane (two lower rows) for the three RMHD 
calculations: MU2 (left), MU10 (middle), and MU200 (right). Scales are logarithmic, and column density N is given in cm 2 . 

These differences in the temperature distribution are fundamen- tude (e.g., iLenzuni et al.lll995t) . The high temperature region is 

tal to the SED estimates. We ran the MU2 calculations until time then able to radiate efficiently and the temperature increases fast 

to + 1.2 kyr, where T max ~ 1 160 K and the outflow extends to towards 2000 K, at which time H2 dissociation starts. 
~ 850 AU, i.e., just before the start of the second collapse. At 
temperature T > 1200 K, dust begins to evaporate, which con- 
siderably lowers the opacity by more than two orders of magni- 
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Fig. 3. Spectral energy distributions for the three models at different times: MU2 (top row), MU10 (middle row), and MU200 
(bottom row). The first panels in the left column show the SEDs obtained from the initial conditions. SEDs are computed at each 
epoch for three different inclination angles: = 0° (pole-on view, solid line), 6 = 60° (dotted line), and 9 = 90° (edge-on view, 
dashed line). 



3.2.2. The intermediately magnetized case - yu = 10 



3.2.3. The very weakly magnetized case - /j. = 200 



In the intermediately magnetic field strength model, there is 
no fragmentation, but an outflow of extent ~ 300 AU and veloc- 
ity ~ 1 km s~' is launched and a rotationally supported disk of 
radius ~ 1 00 AU is formed. This r esult a grees well with previous 
studie s of iHennebelle & Ciardl {2009) and ICommercon et al.1 
d2010l) . At time t + 1.81 kyr (with t = 35.7 kyr), there are 
1.3 x 10 6 cells in the whole computational domain. The mass of 
the FHSC is ~ 0.06 M and its mean temperature is ~ 490 K 
(Tmax ~ 810 K). From the temperature maps, we distinguish two 
shocked regions (i.e., where there is a weak spike in gas tem- 
perature) in addition to the accretion shock on the central FHSC 
itself. In the xy-plane, the shocked region corresponds to the bor- 
der of the disk. In this optically thin region, it is then easy to see 
the spike in gas temperature that is characteristic of radiative 
shocks and has a length o f a few photon mean-free paths (e.g., 
iMihalas & Mihalasll 19841) . In the xj-plane, we identify a shock 
at the outflow border, with a spike in the gas temperature of the 
same amplitude as at the disk border. In comparison with the 
MU2 model, the outflow is broader and denser, and the vertical 
extent of the pseudo-disk is strongly diminished. The calcula- 
tions of the MU10 model were run until T mlLX ~ 1200 K at time 
to + 2.94 kyr, just before the second collapse is expected to begin. 



In the very weakly magnetized model, the collapsing core 
forms a prestellar disk that fragments into several objects (five 
at time to + 3.26 kyr), but no outflow is launched because of 
the weak initial magnetic field. The disk radius extends up to 
~ 150 AU. At time t + 3.26 kyr (with t = 35.4 kyr), there are 
1 .95 x 10 6 cells in the whole computational domain. The mass of 
the fragments is about the same: M ~ 0.043 M Q for the central 
fragment and M ~ 0.036 - 0.039 M for the orbital fragments. 
The mean temperature is T ~ 270 K for the central fragment 
and T ~ 245 - 263 K for the orbital fragments. Heating of 
the gas around the central fragment is much less efficient than 
in the two previous simulations. In addition, we see from the 
temperature maps in the xz-plane that the radiative shock at the 
vertical borders of the disk is sharper and is located closer to 
the center compared to the MU10 run. The infall velocity and 
thus the shock strength are also higher. Before being accreted 
by a fragment, gas has already lost a part of its gravitational 
energy through that shock, and the accretion shock on the 
fragment is then weaker. The MU200 calculations were halted 
before the onset of the second collapse, the first core lifetime 
being r elatively long for the hydrodynamical case dTomida et al.l 
l2010ah . 



6 



B. Commercon et al.: Synthetic observations of first hydrostatic cores in collapsing low-mass dense cores. 



From this qualitative description, we clearly see that the 
three models are different and relevant to the various conditions 
that may be found at the FHSC stage. In the following, we will 
see that these differences are also expressed in synthetic SEDs 
of these collapsing objects. 

4. Synthetic observations 

4.1. Definitions of classical physical stage indicators 

The SEDs of the models were used to calculate classical ob- 
servational signatures. We monitored for each model the tempo- 
ral evolution of the bolometric luminosity Lboi, the bolometric 
temperature 7boi, and the ratio of bolometric to submillimeter 
luminosity L\, \/L smm . The bolometric luminosity Lboi is defined 



Lboi 



f 

Jo 



And l S v dv, 



(1) 



where d is the distance to the object, S v is the flux density, and v 
is the f requency. The bolome tric temperature T\, Q \ is defined fol- 
lowing lMvers & Laddl (1 19931) . i.e., as the temperature of a black- 
body with the same mean frequency as the observed SED, which 
can be approximated by 



C° vS v dv 
T bol = 1.25 x lO- 11 ^ K. 



(2) 



j °°S v dv 

Following lAndre et"ai1 d 19931) . we also calculated the ratio 



Lbo\/L smm , with the submillimeter luminosity estimated as 



Lsn 



I 



v=c/350fim 



4nd 2 S> 



dv. 



(3) 



These three quantities are commonly used to indicate the phys- 
ical stage of observed sources. Lboi and L\, \/L smm are both used 
to define class boundaries. For our pur pose, only the boun dary 
between Class and Class I matters. [Andre et all (1 19931) de- 
fined Class sou rces as objects with L D oi/L smm < 200, while 
IChen et alJ d 19951) defined Class objects with L bol < 70 K. 
Both classifications are problematic for models because one 
has to define the region over which to integrate the luminos- 
ity. Finally, VeLLOs are defined as objects with an internal lu- 
minosity Lj nt < 0.1 L Q , where L mt is the luminosity of the ob- 
j ect in excess of that sup plied by the interstellar radiation field 
ddi Francesco et alJl2007l) . 

4.2. Spectral energy distribution time evolutionary sequence 

In this subsection, we qualitatively describe the SEDs ob- 
tained for the three models. The fluxes of each SED shown be- 
low were computed and integrated over a lobe size of 3000 AU 
x 3000 AU (about half the initial size of the core) and observed 
at a distance of 150 pc (a typical distance to star-forming re- 
gions like Taurus). These choices give an aperture of ~ 20" for 
our synthetic observations, which compares well with current in- 
strument resolutions. Our choice of fixing the same object size 
for all wavelengths is justified because observations often report 
flux over a given area whatever the beam size. In addition, there 
are many different instruments that report star-forming region 
observations and it is impossible to match all beam sizes in a 
single SED. 

Figures [3] shows the SEDs at different epochs for the three 
models. For each epoch, SEDs are presented for three different 
inclination angles: 8 = 0° (pole-on view, solid line), 8 = 60° 
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', 1 i JM I , V , 



1,1, . ■ ■ 

. . L-1 4, 




0=90 
. 6 = 0° 



MU2 MU10 MU200 
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Fig. 4. Radius at which the optical depth 7j = 1 as a function 
of the wavelength for the three models at the same time as in 
Fig. [2] MU2 (black), MU10 (red), and MU200 (blue). The solid 
line shows optical depth with an inclination of 8 = 90°, and the 
dashed line that of 8 = 0°. The dotted line at R = 10 AU roughly 
represents the FHSC's outer boundary. 



(dotted line), and 8 = 90° (edge-on view, dashed line). We also 
computed SEDs for models at an inclination angle of 8 = 45°, 
but these were indistinguishable from those obtained with an in- 
clination of 8 — 0°. As long as f < ?o in all models (see sec. 13.21 
for the definition of fo), SEDs are indistinguishable from the ini- 
tial one and no emission is found at wavelengths shorter than 30 
fim. This behavior is not surprising because the FHSC has not 
yet formed and the gas is still isothermal at T = 1 1 K. 

Surprisingly, the SEDs of the MU2 and MU200 models 
have similar time evolutionary sequences compared to the MU 1 
model, despite their totally different physical conditions. In the 
MU2 and MU200 models, the departure from the initial con- 
ditions starts when t > to and the flux increases at all wave- 
lengths as temperature increases within the FHSC. The peak 
of the SED, however, does not change after the FHSC forma- 
tion. We also find that a significant amount of radiation is emit- 
ted between 20 /im and 100 /mi in the late evolution of the 
MU2 and MU200 models, in contradiction to previous stud- 
ies (e.g., iMasunaga et al.l Il998b lOmukail l2007h . which found 
no observable emission below 30 - 50 ami. This difference can 
be expla i ned by geometrical effects. iMasunaga et al.l d 19981) and 
lOmukail (|20Q7) assumed spherical symmetry, whereas we con- 
sidered a full 3D structure. In our models, the envelope seen 
with inclination angle 8 < 60° is much less dense than that 
of spherical models. In the same manner, the flux we obtain 
with 8 = 90° is strongly reduced in this wavelength range be- 
cause the envelope has reprocessed the radiation emitted by the 
FHSC. Interestingly, the MU10 model does not show any evolu- 
tionary sequence even at the latest stage before the second col- 
lapse, independent of the inclination angle. In the edge-on view 
(8 = 90°), this lack of variation is not surprising, since the MU10 
model combines a disk and a pseudo-disk structure. The radia- 
tion emanating from the FHSC is thus fully reprocessed. But 
even looking pole-on through the outflow, no flux below 30 /mi 
is seen. From Fig. [2] we can see that the outflow in the MU10 
model is much denser than that in the MU2 model; it acts as an 
envelope in which radiation from the FHSC is reprocessed. 
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Figure [4] shows the radius at which the optical depth t,\ = 1 
(the transition from the optically thin to thick regime) as a func- 
tion of the wavelength for the three models at the same time as 
in Fig. |2] and for two inclinations: 9 = 90° (solid line, edge-on) 
and = 0° (dashed line, pole-on). The dotted line at R = 10 AU 
roughly represents the FHSC outer boundary. The sharp jumps 
observed in the wavelength range between ~ 2 pm and ~ 50 
fim are due to jumps in the opacity itself owing to ice or silicate 
for instance (see Fig. [I}- Because the opacity increases when the 
wavelength decreases, the optically thick radius is larger at short 
wavelengths. For 9 = 90°, the optical depth is about the same 
in the three models, since the pseudo-disk or the disk play the 
same role in obscuring radiation from the FHSCs. The central 
FHSC boundary can thus only be seen in an edge-on view at 
wavelengths > 5 mm. The disk or pseudo-disk features alone 
(typical size of about 100 AU) are nevertheless observable at 
shorter millimeter wavelengths, although these features are not 
reliable diagnostics of the presence of FHSCs, since they will 
not disappear as the second collapse proceeds. For 9 = 0°, the 
radius at which the optical depth reaches unity corresponds to the 
FHSC radius for wavelengths > 100 pm, and even marginally for 
shorter wavelengths 3 pm< A < 100 pm for the MU2 model. At 
inclination angle 9 < 60° (since the = 0° SED closely resem- 
bles the 9 = 60° one), the flux received at wavelengths > 100 
pm traces the FHSC outer boundary, and the accretion luminos- 
ity can escape. This result clearly shows the importance of prop- 
erly resolving the accretion shock on the FHSC with a full RHD 
solver and can explain the differences found in the past litera- 
ture. It also has some important implications for observability 
since an FHSC should then be observable as a point source at 
mid- to far-infrared wavelengths if the inclination is favorable. 

Last but not least, there is a clear FHSC lifetime depen- 
dance on the initial magnetization of the cloud. Whereas the 
FHSC has a lifetime > 4 kyr in the MU200 model, it is shorter 
in the other models as magnetization increases, i.e., ~ 3 kyr 
for MU10 and ~ 1.2 kyr for MU2. The evolution is faster as 
magnetization increases, since the accretion rate increases with 
stronger magnetic braking that t r ansports angular momentum 
outward dCommercon et all 120101 : iJoos et ai]|2012l) . As already 
mentioned, iTomida et alj d2010al) reached a similar conclusion 
about the FHSC lifetime using RHD models of different initial 
mass rotating dense cores. 

4.3. Lboi, Thai, and L\, i/L smm 

Figure [5] shows the evolution of the internal luminosity, i.e., 
the difference between the bolometric luminosity at time t and 
the initial bolometric luminosity (top panels) and the bolomet- 
ric temperature (bottom panels) as a function of time for the 
three models and three inclination angles: 9 = 0°, 60°, and 90°. 
Time to is indicated by the vertical dashed line. The dotted line in 
the bolometric luminosity plots indicates the VeLLO boundary 
ddi Francesco et al1l2007l) . L mt < 0.1 L Q . The bolometric lumi- 
nosity and temperature globally increase after the formation of 
the FHSC for inclination angles 9 < 60°. Only the edge-on view 
gives rise to decreasing bolometric luminosity because of the 
disk or pseudo-disk structure. More interestingly, the luminos- 
ity in the 9 < 60° cases increases relatively slowly and becomes 
greater than 0.1 L Q only at late times, long after the FHSC for- 
mation. FHSCs will thus appear as VeLLOs for most or all of 
their lifetime. 

In the bolometric temperature plots, th e horizonta l dotte d 
line shows the Class 0/Class I boundary of IChen et al. I (fl99l . 
i.e., 7boi = 70 K. In terms of 7boi, FHSCs are indistinguishable 



from starless cores and remain at the lower end of the Class 
7boi range. The same behavior is obtained for the ratio Lb i/£ S mm- 
In summary, these classical definitions cannot help in distin- 
guishing between an FHSC and a more evolved Class source, 
or even between an FHSC and a starless core. 

4.4. FHSC signature 

We have shown in the previous section that classical obser- 
vational signatures are not reliable evolutionary indicators dur- 
ing the early collapse stage. Thanks to the increasing spectral 
coverage and sensitivity of new instruments, however, SEDs are 
increasingly more accurate, and new indicators can be defined 
to distinguish between a collapsing core with an FHSC and a 
prestellar (starless) core. 

Figure [6] shows the SED evolution as a function of time and 
inclination for the MU2 model at stages before (black lines) and 
after (red lines) FHSC formation. The blue cross represents the 
sensitivit y at 24 yum (0.83 m Jy at 3cr/24 s) for the c2d Legacy 
program (lEvans et al.ll2003l) and the horizontal blue line repre- 
sents the inflight sensitivity of the Herschel PACS instrument at 
70 /urn, derived from science observations (~ 4.4 mJy in point 
source mode, 5cr/lh, iPoglitsch etaDl2010t) . Shortly after FHSC 
formation, the SED does not significantly depart from the ini- 
tial one. After the FHSC begins to build up its mass, however, 
the accretion shock becomes stronger and departures from the 
initial SED are much clearer. To take into account the uncertain- 
ties on the initial temperature of the prestellar cores, we show 
SEDs of prestellar cores with initial temperatures ranging from 
9 K to 13 K (about 20 % upper limit uncertainty; Launhardt et 
al., in prep). The SEDs with accreting FHSC are clearly distinct 
from the initial uncertainties on the temperature for inclination 
9 < 60° although the SED peak itself does not show evidence 
of an FHSC. According to the sensitivity of the Herschel PACS 
instrument, FHSCs can be identified with sufficiently deep ob- 
servations. More interestingly, the sensitivity of Spitzer at 24 pm 
could also be sufficient for FHSC detection at late stages of their 
lifetime. The wavelength range over which the flux has increased 
the most after the FHSC formation is 20 pm < A < 100 pm (at 
least a factor 10 at 100 pm for 9 < 60°), whereas L smm increases 
only by a factor less than 2. Similarly, the ratio L\, \/L smm in- 
creases at maximum by only a factor ~ 2 after FHSC formation. 

The estimate of L\, \ from observations strongly depends on 
the number of points in the SED, i.e., the spectral coverage. 
Instead of the physical stage indicator Lboi/^smm, we chose here 
to express a related luminosity ratio using a particular wave- 
length. In this example, we chose the luminosity at 70 pm to 
show how it is possible with the Spitzer/MlPS or Herschel/PACS 
instruments to identify FHSC candidates in comparison to star- 
less cores. We note that a similar analysis can be performed us- 
ing the 100 pm, depending on the available observations. The 
flux at this wavelength is well constrained, and we just showed 
that L smm increases only by a factor 2 after the formation of the 
FHSC. Finally, we note that similar flux increases in the wave- 
length range 20 pm < A < 100 pm are also observed in 10 M Q 
dense core collapse calculations since, even if the envelope is 
thicker, the accretion shock on the FHSC is also stronger (e. g., 
for massive dense core models, see ICommercon et al.|l201 lbh . 

Figure|7](top) shows the time evolution of the flux at 24 pm 
(solid line), 160 /urn (dotted line), and 1200 pm (dashed line) as a 
function of the flux at 70 pm and of the inclination angle for the 
three models. The horizontal gray line indicates the aforemen- 
tioned 24 pm sensitivity limit of Spitzer and the vertical dotted 
line the 70 pm sensitivity limit of the Herschel PACS instru- 
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Fig. 5. Bolometric luminosity (top) and bolometric temperature (bottom) as a function of time for the three models: MU2 (left), 
MU10 (middle), and MU200 (right). Bolometric quantities are computed from models with three different inclination angles: 9 = 0° 
(black), 9 = 60° (red), and 9 = 90° (blue). 
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Fig. 6. Spectral energy distribution evolution as a function of time and inclination for the MU2 model. The black line represents the 
SED when the FHSC has not yet been formed, whereas the red lines indicate SEDs after FHSC formation. The gray area around the 
initial SED reflects the uncertainties in the initial temperature of the prestellar c ore (9 K< Tp <13 K ). The blue cross represents the 
sensitivity at 24 fim (0.83 mJy at 3cr/24 s) for the Spitzer c2d Legacy program dEvans et alj|2003l) and the h orizontal blue line the 
sensitivity of the Herschel PACS instrument (~ 4.4 mJy in point source mode, 5cr/lh; Poglitsc h et al.ll2.Q10h . 



ment. For < 60°, the flux Fjq increases by about two orders 
of magnitude after FHSC formation for the MU2 and MU200 
models, whereas the MU10 model does not show any increase, 
as expected, since its SED does not show any evolution. Note 
that the flux received at 70 //m is always above the sensitivity 



limits of the current instruments. The most dramatic increase 
is found in the flux at 24 fim, which is almost zero at FHSC 
formation and then increases to detectable values at late times. 
Similar to L smm , the fluxes at 160 fim and 1200 fim remain al- 
most constant throughout the FHSC lifetime. The fluxes at 24 
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Fig. 7. Top: Evolutionary diagram of the flux at at 24 pm (solid line), 160 fim (dotted line), and 250 fim (dashed line) as a function of 
the flux at 70 fim for the three models and for the the same three inclination angles as previously. The horizontal gray line indicates 
the aforementioned 24 pm sensitivity limit of Spitzer and the vertical dotted line the 70 pm sensitivity limit of the Herschel PACS 
instrument. Bottom: Color-color plot representing log(/lF / | 1 60/ AF,{10) as a function of log(/lF / |70/ XP^2A) for the three models with 
the same three inclination angles. 



fim and 70 fim can thus help in identifying FHSC candidates 
with respect to starless cores if the physical conditions are ap- 
propriate (either a strong or a very weak magnetic field), and 
if the objects are observed with a favorable inclination < 60°. 
The evolutionary tracks of the FHSCs in the color-color plane of 
log(/lF /) 70/^24)-log(/lF / ) 160/^70) are shown in Fig.|7](bot- 
tom) for the three models at the three same inclination angles 
as discussed previously. There is a clear evolutionary track fol- 
lowed by the accreting FHSC that can be directly compared to 
observations. 



5. Discussion 

5. 1 . Comparison to previous work 

Spectral energy distributions of FHSCs embedded in their 
parent cores have already been studied for instance b y 
Boss & Yorkel dl995h . IMasunaga et aD dl998h. lOmukail (I2007I) . 

Saigo & Tomisakal d20l¥ . While 



Tomida et alJ d2010ah . and 



they all found that the radiation emitted by the FHSC is es- 
sentially completely reprocessed by the dust in the surround- 
ing core to 50 - 200 fim with no observable emission be- 
low 30 fim, we found that there is indeed observable emis- 
sion at shorter far-infrared wavelengths down to 20 fim. These 



differences may be explained by the different physical mod- 
els and post-processing of each study. W hile we directly post- 
processed 3D AMR RMHD calculations. iBoss & Yorkel dl995h 
used a spherical grid and integrated the RHD equations with a 
low numerical resolution. Their FHSC temperature is also rela- 
tively low (< 200 K). Their results are therefore similar to ours 
shortl y after FHSC formation . In a ddition, a s men tioned previ- 
ously, IMasunaga et alJ d!998l) and lOmukail d2007l) used spher- 
ical model s and their SEDs a r e then similar to ours with 6 = 
90°. Next, ISaigo & Tomisakal d2011l) post-processed 3D high- 
resolution collapse calculations with plane-parallel approxima- 
tion but neglected magnetic fields and used a barotropic EOS. 
They observed an increase of the total energy in the SED but 
their SEDs were calculated over a very small area of 40 AU x 
40 AU and therefore excluded t he large-scale emission dominant 
at wavelengths > 100 fim. Last. lTomida et al.l d2010al) computed 
SEDs from 3D RHD calculations but their cores do not frag- 
ment (slow rotation). In addition, they used Bonnor-Ebert initial 
density profiles in which accretion rates are lower and therefore 
FHSCs are fainter. 
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Fig. 8. Dust continuum emission maps for the MU2 models at four wavelengths, from left to right: 24 fim, 70 fim, 160 fim, and 
1200 fim. Three different times are shown that trace the evolutionary stage of the source: to - 26.43 kyr (top row, i.e., starless core), 
fo + 0.7 kyr (middle row, i.e., faint FHSC) and to + 1.13 kyr (bottom row, i.e., end of the FHSC lifetime). The fluxes are given in 
Jy/P x (1 P x = 0.08" x 0.08") for objects at a distance of 150 pc and observed pole-on. 



5.2. What can SEDs tell about physical processes and 
evolutionary stages? 

Figures|3]and|6]show that the evolution of an FHSC through- 
out its lifetime can be reflected in observable changes of the 
SED, but only under certain conditions. Only the intermedi- 
ate MU10 model does not show any evolutionary sequence be- 
cause of the thick envelope structure consisting of a system 
disk + pseudo-disk + outflow. The MU2 and MU200 models 
show an evolutionary sequence very different than that of the 
MU10 model, but very similar to each other. The MU2 model is 
strongly magnetized and exhibits a pseudo-disk + outflow sys- 
tem, whereas the MU200 model fragments and exhibits a ro- 
tationally supported disk. An FHSC with a clear evolutionary 
sequence can then be observed in simple geometrical config- 
urations and with inclination angle < 60°. Nevertheless, it 



is impossible to distinguish between a magnetized and a non- 
magnetized scenario from a strict SED point of view. High an- 
gular resolution images of, at least, dust emission are required to 
distinguish between these scenarios. This kind of diagnostic can 
only be provided by interferometry at longer wavelengths (see 
Paper II). 

We have shown than it can be relatively straightforward to 
distinguish a starless core from an FHSC thanks to flux re- 
ceived at mid- to far-infrared wavelengths, i.e., at 24 fim and 
70 fim, but it remains hard to distinguish be tween an FHSC an d 
a young stellar object (YSO) with a disk dPineda et al.ll2.Ql il) . 
i.e., an accreting second hydrostatic core (SHSC). This difficulty 
is due to the thick envelope in which the FHSC and SHSC are 
each embedded throughout the Class phase that reprocesses 
short-wavelength radiation to longer wavelengths. Here, mod- 
els clearly fail to predict any detectable differences. Using only 
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dust continuum emission is not the appropriate way to differen- 
tiate the two core stages, since almost all radiation emanating 
from both the FHSC and the SHSC will be reprocessed by the 
thick envelope. Another approach could consist of focusing on 
the outflow properties (collimation, velocity, momentum, etc.) 
and on the chemical properties of the sources, such as gas phase 
and ice mole cules abundance s, which reflect the local thermal 
history (e.g., iKim et aU201 lb . We will study these properties in 
upcoming papers. 

5.3. Example of methodology to target first core 



We showed in Sec. 14.41 that the fluxes at 24 fim and 70 fim 
can be good indicators of FHSCs in comparison with starless 
cores. We also did not find any significant emission below 10 
fim. We infer that the combination of a point source detection at 
70 fim or even at 24 fim with no detection at wavelength < 10 
fim can be a good indicator of FHSCs from a strict SED point 
of view. Under the current telescope sensitivity limits, a detec- 
tion at 70 fim alone is most likely, since the emission at 24 fim is 
faint, and reaches the Spitzer sensitivity limit only in the MU2 
model in the second half of the FHSC lifetime. Figure [8] shows 
dust continuum emission maps for the MU2 model at three dif- 
ferent evolutionary stages and for four wavelengths (24 fim, 70 
fim, 160 fim, and 1200 fim). It clearly shows that the emission 
at 70 fim and 24 fim should be point-like after the FHSC for- 
mation. Note that at all evolutionary stages, a flux at 70 fim can 
be detected but its emission would not look like a point source 
for a starless core. It will instead be much more extended (i.e., 
the size of the initial dense core) because of the ISRF repro- 
cessing. Thanks to the recent Herschel and Spitzer observations, 
observers should be able to combine data from both telescopes 
and make comp arisons with our predicted color-color plot in Fig. 
[7] For instance. Idi Francesco et al.l d2010l) already used this pro- 
cedure and combined maps of Herschel data with Spitzer ones 
to distinguish between prestellar and protostellar cores. Once 
FHSC candidates are identified from these flux examinations, 
the next step will be to target them with high-resolution interfer- 
ometry such as ALMA (see Paper II) to constrain the physical 
properties (e.g., strength of the magnetic field). 

As mentioned earlier, a complementary approach would be 
to use chemical analysis to separate the FHSC and the SHSC 
stages. At temperatures above ~ 1000 K, grain evaporation starts 
and refractory chemical species such as C and Si are released 
into the gas phase, cons iderably changing th e chemical network 
and the opacities (e.g., iLenzuni et al.ll 1 995b . Deuterium-related 
chemistry could also trace an SHSC because deuterium-bearing 
species should be unde r-abundant in high-temperature regions 
dAlbertsson et alj|2011l) . 

Another more direct and more promising way to stage dif- 
ferentiation would be to find evidence of a low-velocity out- 
flow without its high-velocity counterparts. Recent theoreti- 
cal studies have shown that a low-v elocity outflow can be 



launched during the FHSC stage (e.g., Hennebelle & Fromanj 
120081: ICommercon et al.ll2010HTomida et al.ll2010bl) without the 
high-velocity outflow that is always present at the SHSC stage 
dMachida et al.ll2008l) . Not only may the morphology and ener- 
getics of the outflow help in identifying the cores, but also the 
chemical composition of the gas in the outflow. Since tempera- 
tures are different during the FHSC and SHSC stages, we may 
expect that the low-velocity outflow would show evidence of 
low-temperature gas-grain chemistry and that the high-velocity 
outflow would show evidence of high-temperature chemistry. In 
addition, the low-velocity outflow is launched from regions with 



temperatures < 1000 K where grains are still not fully evapo- 
rated. Grains should then be present in the FHSC outflow and 
since the velocity is relatively low (a few km s~'), they are not 
expected to be destroyed in the shock at the envelope interface. 
Therefore, we may expect FHSC outflows to be seen in thermal 
dust emission if they are dense enough. It is worth mentioning 
that some dust continuum observations may be also contami- 
nated by CO emission from o utflow gas, making them easier to 

see (e.g., iDrabek et al1l2012h. 

Recent work bv lPrice et al.l d2012l) reported the formation of 
a collimated jet driven from the FHSC, support ing recent ob- 
serva tions of an FHSC candidate in Per-Bolo 58 dDunham et alj 
1201 lb . The timescale needed to produce such a collimated jet (~ 
2 kyr), however, is longer than the expected lifetim e of FHSCs 
in a m agnetized collapsing dense core. In addition, iPrice et all 
d2012l) used sink particles (scales smaller than 5 AU are not de- 
scribed) and a barotropic approximation with an adiabatic index 
of y = 7/5, for which FHSC lifetimes are even shorter than those 
we predict in this study with y — 5/3 (the FHSC contracts adi- 
abatically faster wi th y - 7/5 ) . It is therefore very likely that 
FHSCs reported in IPrice et all d2012l) should have already un- 
dergone second collapse and formed an SHSC within ~ 2 kyr 
after their formation. 



5.4. Limits of the model 

In this study, we used the gray radiative transfer approxi- 
mation for the RMHD calculations which may not be the most 
appropriate method at l ater stages of the e volutions. As men- 
tioned earlier, however, IVavtet et"aTI d2012l) showed that multi- 
frequency collapse calc ulations were well - reprodu ced by the 
gray FLD calculations of ICommercon et all (1201 lab , at least for 
the first collapse and FHSC stages. In addition, we used a con- 
stant ratio of specific heats y, set to 5/3, although it is known 
that y should change to a value of 7/5 as temperature increases, 
and H2 rotational d egrees of freedom become available (e.g., 
iMachida et"aDl2008h . With y = 7/5, FHSCs compress and cool 
more efficiently and thus reach the second collapse more quickly 
(e.g., ICommercon et al.ll2.010h , which reduces the FHSC life- 
time. 

For the RAD MC-3D post-processing , we used one specific 
opacity model of lSemenov et al.l d2003l) for all temperatures and 
densities. This choice may not be appropriate since grain man- 
tle compositions change as temperature increases. Nevertheless, 
the changes in opacity remain relatively minor at temperatures 
< 1000 K and we do not expect them to have a significant impact 
on our results. Finally, we have neglected scattering although 
it has been found to be important for wav elengths < 10 fim 
dYoung & Evansll2005l: iDunham et ail 1201 Ol) . Since we did not 
find any emission from FHSCs at these wavelengths and we do 
not consider this wavelength range significant for our conclu- 
sions, scattering effects are not relevant for our purposes. 



6. Summary and perspectives 

We used state-of-the-art radiation-magneto-hydrodynamic 
models of the first collapse and first hydrostatic core formation 
in the context of low-mass star formation to predict spectral en- 
ergy distributions and related observational indicators. We in- 
vestigated the collapse and the fragmentation of a 1 M dense 
core with three different magnetization levels up to the start of 
the second collapse. We showed that the FHSC lifetime strongly 
depends on the magnetization of the collapsing core and that a 
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strong magnetic field leads to a shorter lifetime because mag- 
netic braking raises the accretion rate. We post-processed the 
RMHD calculations with a 3D radiative transfer code and pro- 
duced synthetic spectral energy distributions. We showed that 
FHSCs can be identified thanks to the rapid increase in the 
flux received between 20 /mi and 100 //m for inclination an- 
gle 9 < 60° under different physical conditions (e.g., high or 
weak magnetization). We also found that FHSCs are globally 
consistent with VeLLOs. The classical stage evolutionary indi- 
cators 7boi, Lbo\, and Lboi/^smm are not reliable for identifying 
FHSCs since they do not show any significant evolution through 
the FHSC lifetime. We found that a point source detection at 
wavelengths shorter than 200 /mi is a good indicator of FHSC in 
combination with the absence of detection at wavelengths < 10 
/zm . We currently cannot, however, provide any diagnostics to 
distinguish an FHSC from a low-luminosity SHSC. 

Recently, Herschel has mapped hundreds of dense cores 
through, e.g., the EP OS and Gould Belt Survey key observa- 
tion programs (e.g., iHenning et all 120101: lAndre et al.l 1201 Ot 
iBontemps et al.l |2010|) with PACS data at 70 //m. The com- 
munity has also gathered many Spitzer observations at 24 //m 
and 70 /mi. Combining the observations for common sources 
of both datasets could identify FHSC candidates. Cataloging 
these candidates will provide the basis required for observation 
proposals aimed at high-resolution interferometers such as the 
Submillimeter Array (SMA) and ALMA, which will be indis- 
pensable to accurately measure the physical properties of the 
collapsing cores (see below). Future progress is also expected 
from the next generation of space telescopes such as SPICA with 
the SAFARI instrument (with 1-sigma r ms sensitivities of ~ 2 
/Jy at 24 //m and ~ 50 //Jy at 70 jum; ISwinvard et~ai]l2009t) . 
which could then test the evolutionary tracks we predicted. 

SED information alone cannot help in distinguishing the im- 
portant physical mechanisms at work at early Class stages 
(large disk formation and fragmentation, outflow launching, 
etc.). Only interferometric observations would be capable of 
such a feat, and forthcoming ALMA observations will give a 
clearer, if not definitive, answer to the fragmentation properties 
of Class sources. ALMA dust continuum emission synthetic 
maps will be presented in Paper II. Predictions other than dust 
continuum emission only are also needed to target FHSC candi- 
dates more efficiently. Chemistry post-processing and molecular 
line emission calculations are therefore a natural next step for 
FHSC observation predictions. Future investigations will also 
aim to follow the second collapse with the appropriate physics 
(mainly H2 dissociation). 
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MU10, Resolution 
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Fig. A.l. Spectral energy distribution for the MU 10 models com- 
puted on three different grids: the AMR grid for the black lines, a 
uniform 256 3 cube for the red lines, and a uniform 5 12 3 cube for 
the blue lines. Two different inclination angles are also shown: 
6 = 0° (solid line) and 9 = 90° (dashed line, almost indistin- 
guishable from one to another). 



Appendix A: AMR versus cartesian grid for 
postprocessing 

In this appendix, we present SEDs computed with two dif- 
ferent post-processing methods. The first is the one we used 
throughout this paper, i.e., the AMR grid of RAMSES is directly 
loaded in RADMC-3D. The second projects the AMR data of 
RAMSES into a uniform cube of extent ~ 3000 AU x 3000 AU. 
The second method is much simpler for interfacing, but a lot of 
information is lost when smoothing the finer levels onto the uni- 
form grid. In our calculations, since the finer levels correspond 
to the high-temperature and high-density regions, smoothing to 
a lower resolution may affect the SEDs. 

Figure [AT] shows SEDs for the MU10 models computed on 
three different grids: the AMR grid for the black lines, a uniform 
256 3 cube for the red lines, and a uniform 5 12 3 cube for the blue 
lines. Two different inclination angles are shown: = 0° (solid 
line) and 9 = 90° (dashed line). The three SEDs obtained with 
9 = 90° are similar with no emission below 30 fim independent 
of the grid used for post-processing, but they are different for 
9 = 0°. For the latter case, there is indeed emission between 3 fim 
and 30 fim and an increase of the flux around 100 /mi. We clearly 
see that poor resolution may strongly affect the SEDs, which 
solely comes from the fact that high-density regions have been 
smoothed out by reprojection of the AMR grid onto a uniform 
one. 



